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QUANTUM GALVANOMAGNETIC AND THERMOMAGNETIC EFFECTS IN 
GRAPHITE TO 18 TESLAS (180 kG) AT LOW TEMPERATURES 
by John A. Woollam 
Lewis Research Center 

SUMMARY 

The electrical magnetoresistance, Hall resistivity, adiabatic thermoelectric power, 
adiabatic Nernst-Ettings hausen coefficient, and thermal resistivity are measured as a 
function of magnetic field from 0 to 18. 3 tesla (1 T = 10 kG). Temperatures range from 
1. 1 to 4.2 K. From the measured coefficients, the isothermal thermopower, electrical 
conductivity, and Hall conductivity are found. From a study of these effects, an assign- 
ment of electrons and holes to positions in the Brillouin zone is made. Spin splitting of 
the n = 1 electron Landau level is found in three of the observed coefficients. The 
thermopower quantum oscillations at low temperature and high magnetic field have dis- 
torted shapes. These shapes are compared with predictions of various theories, and 
are found to be correctly predicted by a theory by Horton. Sugihara and Ono have pre- 
dicted positions in field for Landau level crossings, and we find good agreement with 
their predictions with the exception of the n = 1 electron level. A study of deHaas - 
van Alphen periods and anisotropy shows that majority carriers are the same as found 
in natural single crystals. 


INTRODUCTION 

Pressure annealed pyrolytic graphite is a highly ordered form of carbon (ref. 1) but 
is not a single crystal. The lattice for a single crystal is hexagonal. For pyrolytic 
samples there are individual crystallites with [0001] axes nearly parallel to each other, 
but [1100] axes, for example, are randomly oriented to each other. Typical crystallite 
size is 10- to 100-micrometer diameter. Because of the nature of the Fermi surface of 
graphite, the properties of pyrolytic and natural single crystals are very similar, as 
shown in this and other papers (ref. 2). 

Fermi surface studies of natural single crystals of graphite were made by Soule, 


McClure, and Smith (ref. 3) in fields to 2.4 tesla (1 T = 10 kG) using the Shubnikov- 
deHaas effect. deHaas ~ van Alphen studies on pressure annealed pyrolytic graphite were 
made by Williamson, Foner, and Dresselhaus (ref. 2). Magnetoreflection results 
(ref. 4) were used along with the deHaas - van Alphen effect to determine the magnitude 
of the parameters necessary to describe the energy band structure derived by 
Slonczewski and Weiss (SW) (ref. 5). 

Soule studied the effect of boron doping on the deHaas - van Alphen frequencies 
(ref. 6). Boron adds positive charge and thus shifts the Fermi energy and changes the 
cross-sectional areas. The shifts in majority carrier frequencies assign the location 
of electron and hole carriers in the Brillouin zone and thus determine the signs on sev- 
eral of the parameters in the SW (ref. 5) band theory. Recently Schroeder, Dresselhaus, 
and Javan (ref. 7) have used magnetoreflection to determine the band parameter signs 
and find a conflict with the assignment found by majority carrier deHaas - van Alphen 
frequency shifts. 

In this investigation two independent methods were used to determine these signs: 
Quantum extrema in the Hall coefficient are compared with a theory by Argyres (ref. 8) 
and from this a direct determination of electron and hole carrier positions in the 
Brillouin zone is made. This then determines the signs of the band parameters. Quan- 
tum resonances in the thermopower are used as a second and independent method of de- 
termining carrier locations and band parameter signs. 

Wagoner (ref. 9) studied the electron spin resonance as a function of temperature and 
determined an effective g factor, which is an average over all carriers in the Brillouin 
zone for particular magnetic field directions. We have made the first observation of 
spin splitting of Landau levels in graphite. These determine the splitting at particular 
points in the Brillouin zone rather than an averaged value. 

In this report the first studies of quantum thermomagnetic effects in graphite are 
made, as well as Shubnikov-deHaas studies in pyrolytic graphite. Temperatures range 
from 1. 0 to 4. 2 K and fields up to 18. 3 tesla. Thermopower quantum oscillations have 
highly distorted line shapes and these are compared with various theories. The 
Shubnikov-deHaas results are compared with the theories of Adams and Holstein 
(ref. 10). Sugihara and Ono (ref. 11) have calculated the positions in field for Landau 
level crossings of the Fermi energy and show that at high fields the crossings are not 
predicted by strict periodicity in inverse magnetic field. Experimental and theoretical 
crossings are compared for each of the observed effects. 


THEORY 


Kinetic Equations and Transport Coefficients 

The linear relations between current density J, electric field E, negative temper- 
ature gradient G, and heat current density w, may be written (ref. 12) as 




E* =p’J + e’w* 

Adiabatic 

(i) 



G = n'J + yw* 


(2) 

or independently as 






E* = p J + eG 

Isothermal 

(3) 



w* = — 77 J + Xg 


(4) 

where 


E* = E - 

e 


(5) 



W* = W - 

e 


(6) 

and /± 

is the chemical potential. 

From equations (1) and (3) when J is zero 




e’w* = eG 


(7) 

where 

e is the isothermal thermoelectric tensor, and e 1 

is the adiabatic thermoelec- 


trie power. 

Equation (2) shows that when J = 0 


G = yw* (8) 

so 

e’ = ey (9) 
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where y is the thermal-resistivity tensor and 


T— l 
1 

<✓< 

II 

(10) 

A A ^ A 

e = e'A 

(ID 


where A is the thermal- conductivity tensor. 

In the present experiments transport coefficients were measured as a function of 
magnetic field strength, temperature, and angle between the magnetic field and the C 
axis of graphite. The specific coefficients measured were the yy and yx components 
of the p, e’, and y tensors. In addition, the following coefficients were calculated: 

°yy and CT xy fr ° m 

a = p -1 (12) 


A X y and Ayy from equation (10), and from equation (11). It is important to point 

out that the coefficient known as the isothermal thermopower, is the negative of the 

absolute thermopower or thermopower S, defined by some authors (e. g. , Ziman; 
ref. 13). In this report, S is thermopower, not e^. 


Adams and Holstein Theory for o yy 

Adams and Holstein (ref. 10) calculate the electrical conductivity for the case of 
crossed electric and magnetic fields for both degenerate and nondegenerate solids. Per- 
tinent to this study is their calculation of the conductivity for degenerate systems in the 
limit of large magnetic fields and low temperatures. Their calculation is quantum me- 
chanical in that a density matrix formalism is used. 

The calculation considers various types of scattering mechanisms. At zero degrees 
the conductivity, valid for any particular scattering mechanism, is characterized by 
singularities as a function of magnetic field: 


a 


yy 






Terms in the square of 
_the scattering potential. 


(13) 


where n and m are Landau level quantum numbers and is the Fermi energy. The 
cyclotron frequency w is given by 


w = 


eH 

m*c 


(14) 
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where H is magnetic field, and m* is the effective mass. The term in brackets is de- 
termined by the type of scattering. Equation (13) has infinite discontinuities where the 
Fermi energy and the Landau level coincide, that is, when 



At nonzero temperatures, there is thermal and collision broadening of the Landau levels 
and equation (13) leads to an oscillatory part, a , of the conductivity given by 

(2 7 rM)- 1 / 2 cos(2EZ-l- 7 rM) (16) 

V H 4 / 

where is the nonoscillatory part of and F is the deHaas-van Alphen fre- 

quency, proportional to the extremal cross-sectional area of the Fermi surface perpen- 
dicular to the magnetic field. The frequency F is related to S p and cj by 

F _ g F 
H hco 

Equation (16) shows that the conductivity oscillates periodically in 1/H. 

In graphite, the effective masses are so small (ref. 3) that for fields above 1 tesla 
terms other than the M = 1 term in equation (16) become important, and the oscilla- 
tions become highly distorted. Another important feature of equation (16) is that for low 
temperature and high magnetic field the positions of the extrema in a will shift. At 
T = 0 K, a discontinuities occur at coincidence of the Fermi energy gp and the 
Landau level. At higher temperatures the extrema move to lower fields. It is also im- 
portant to know that a does not involve the charge of the transport carrier (ref. 9). 
Thus Oyy is a maximum at coincidence of £p and a Landau level, independent of 
whether carriers are electrons or holes. 



Hall Conductivity Theory 

Argyres (ref. 8) calculates the Hall conductivity cr^ considering only elastic scat- 
tering (ionized impurity scattering, for example) in the Born approximation. The effect 
of both magnetic and electric fields on the collisions is considered, and, as Adams and 
Holstein (ref. 10) have done, he uses the density matrix formalism. For high symmetry 
directions of the magnetic field in graphite (parallel to the C axis, for example) 
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yx 


= -a. 


xy 


where a is the Hall conductivity (ref. 12). Under these conditions Argyres results 

xy 

lead to 



T n[ f n T i 1(n ) + f n + l T i 1 ( n+1) ]| 


2 2 , 
" T n + 1 


(17) 


for isotropic scattering at T = 0 K. In equation (17) q is positive for holes, negative 
for electrons and has the magnitude of electronic charge, and 



(18) 


and 


t" 1 = rjV) + t'V + 1) 


(19) 



( 20 ) 


In these equations f is the Fermi distribution function evaluated at the energy of the 
n n Landau level and f^ is the energy derivative of the distribution function. 

Equation (17) shows that cr xy has positive maxima at coincidence of and the 
Landau level for electron carriers. For hole carriers cr xy has negative extrema at 
coincidence of f F and the Landau level. Thus, Hall effect quantum extrema, as a 
function of field, depend on charge carrier sign while conductivity extrema do not. 

Horton (ref. 14) has used the density matrix approach to calculate the Hall conduc- 
tivity at low temperatures and high magnetic fields. Horton's results are valid for non- 
zero temperatures and predict the same general result as Argyres when temperatures 
are low: cr xy has negative extrema at coincidence of j? f and the Landau level for holes. 


Quantum Thermopower Theory 

In this section we present three approaches to the theory of quantum thermopower, 
in order to compare their applicability to experimental observations. 

Generalized Mott equ ation. - The diffusion thermopower in the presence of a mag- 
netic field is given by (ref. 15) the generalized Mott formula 
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( 21 ) 


3q L d! -L 

fe F 

where S, p, and a are field dependent tensors. As before, the factor q is positive for 
holes, negative for electrons, and has the magnitude of the electronic charge. By defi- 
nition (refs. 12 and 13) 

S(H) = -e(H) (22) 

where e is defined by equation (3). Equation (21) is a result of the linearized Boltzmann 
equation with elastic scattering. In graphite we find experimentally 


and 


l°xyl << l a yyl 


From symmetry, 


Thus, 


0(T 


3ct , 

xy 

« 

yy 

3H 


3H 

= a 

= or 

= cr„ 

yz 

zx 

zy 


s yy (H) 


7T 2 k 2 T 

3q 



(23) 


(24) 


(25) 


From equations (13) and (20) combined with (25) for S , we see that near coincidence 
of i p with a Landau level 

S °c »Vt 6 -l ( 26 ) 

yy 3qjP F 

where 6 - ^ has positive spikes at coincidence of and a Landau level. Thus 
will have positive spikes for hole carriers (+q) and negative spikes for electron carriers. 
For T ~ 2 K and above, the spikes in graphite become rounded maxima. Close to 4 K 
in graphite, is purely sinusoidal in 1/H. 

Thermodynamic theory for S yy . - To approximate the amplitude and line shape of 

the spike shaped oscillations that occur at low temperature and high field, we tried an 
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argument discussed by Callen (ref. 16), and by MacDonald (ref. 17). They find that in a 
steady-state condition, and where the Joule heating is small, the Thomson heat is the 
specific heat of the electrons per charge, as a first approximation. This assumes a re- 
versible process or quasi-equilibrium. It then follows that the thermopower, S yy is 



(H) 


a(H) 

Nq 


(27) 


where Q!(H) is the entropy of the system of electrons given by 

a(H) = - MM 
9T 


(28) 


where ft(H) is the free-energy of the system. Equation (27) was also derived by 
Obratzof (ref. 18) and others (refs. 19 and 20), but they assume that 


xy 1 


» k 


yy' 


a condition we did not find in graphite. 

Lifshitz and Kosevich (ref. 21) have evaluated fl under conditions of strong quan- 
tization, that is, for 


hco » kT 

Using the Lifshitz-Kosevich expression for fi, and evaluating equations (27) and (28) we 
find 



2k 

^o 



-MX. 


D 


M 


3/2 


Ag(MX) cos 


/2ttMF 

V H 


27rMy T l N \cos 7rM g m * 

4/ 2m 0 


(29) 


where n Q is the density of particles, A is the extremal cross-sectional area of the 
Fermi surface perpendicular to the field, k is the component of wave vector parallel 
to the field, y is a phase factor normally equal to 1/2, g is a spin factor, and 
is the ’’Dingle” factor which accounts for collision broadening of Landau levels (ref. 22). 
The term A^(MX) is given by 

A, (MX) = 1 ~ MX coth MX (30) 

6 sinh MX 
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I 


where 


27r 2 kT 

Kco 


(31) 


The term Ag for M = 1 is shown for graphite as a function of temperature at various 
fixed magnetic fields in figure 1, and as a function of temperature at 1 tesla for M = 1, 
2, 3, and 4 in figure 2. From figures 1 and 2 and equation (29), it is obvious that the 
factor Ag determines the amplitude and harmonic structure of the oscillations as a 
function of magnetic field and temperature. At low temperatures and high fields the 
M = 1 term in equation (29) decreases in relation to the M = 2, 3, etc. , terms. This 
causes a highly distorted oscillation shape at low temperatures and high fields. A plot 
of equation (29) with parameters applicable for majority carrier holes in graphite at 2 K 
is shown in figure 3. Figure 3 is a plot of Syy as a function of field showing distorted 
oscillation shapes at high fields. Equation (29) predicts better than order of magnitude 
agreement with experimentally observed amplitudes. However, the shape of the dis- 
torted peaks is not correctly predicted and this result will be discussed later. It is im- 
portant to present the thermodynamic, Obratzof, etc. theories because the present ex- 
periments test the applicabilities of these theories to graphite. 


Magnetic 



Figure 1. - A 3 (MX) as function of temperature (eq. (30)) for several field 
strengths. Effective mass, 0.04 appropriate for holes in graphite. 


9 


i. 


Thermopower amplitude determining function, A 3 OVIX) 



Figure 2. - A 3 (MX) as a function of temperature (eq. (30)) 
for several M values, at 1 tesla. Effective mass 0.04 m 0 
appropriate for holes in graphite. 


Thermodynamic thermopower, S W (H), arbitrary scale 



EM' 


M- 


3/2 


A 3 (MX) cos 

L H 



Figure 3. - Theoretical thermodynamic thermopower 
{eq. (29)) as function of field at 2.0 K for hole 
carriers in graphite. 


Horton thermopower theory . - As will be discussed later, the experimentally ob- 
served distortion of oscillation peaks in thermopower is correctly predicted by a theory 
due to Horton (ref. 14). Horton's calculation is quantum mechanical in that he uses the 
density matrix formalism. Scattering is assumed to be due to random point impurities, 
and the calculation is for free electrons. His final expression is 


g _ 2ir 2 n o qT k 2 T 
yy ~ 3 m 


1 - 


15 




8tt 



xM 


A,(MX)sin Mb 

Vm 


where 



A 4 (MX)cos Mb 


(32) 


a 4 (mx) 


— csch(MX) [2 coth(MX) - 2 MX coth 2 MX + MX 
2 L 


and Ag is given by equation (30). The letter q is used rather than Horton's e, n Q is 
the charge carrier density, r the average time between collisions, and 

b = (33 ) 

H 4 

Figure 4 shows S from equation (32) plotted at T = 1.2 K. Figure 4 is plotted 
for q + (hole carriers), and there are higher positive spikes than negative dips. In 
part (a) a Dingle factor (eqs. (29) and (31)) was assumed with Tp = 1.0 K. In part (b) 
Tp = 0. 5 K. As will be seen later, figure 4(b) closely approximates the experimental 
results. Williamson measured a Tp from deHaas - van Alphen effect of 1. 5±0. 1 K for 
a similar sample of pyrolytic graphite. This agreement supports the use of a Dingle 
factor in Horton's expression (introduced in eq. (29)). 

Notice the highly distorted oscillation shapes in figure 4. For the same parameters 
the deHaas - van Alphen effect is not nearly as distorted. There are two reasons for the 
high distortions in thermopower. First, the Ag term shown in figures 1 and 2 causes 
large harmonic content; that is, the inclusion of many of the M terms in equation (32) 
is necessary. The second reason is that there are two terms in equation (32). At low 
fields the sine term is larger, but at higher fields and low temperatures the cosine term 
causes significant distortion. 
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%- t5 »“Z$v{ff VnWH 

n 

+ 79. 5 Hl'2VHge- 27,2nT D ,H A3Sin[f -?] 




(b) Dingle temperature, 0.5K. 

Figure 4. - Theoretical Horton thermopower (eq. (32)) as function of 
field at 1.2 K for hole carriers in graphite. 


The Energy Band Structure of Graphite 

We have found, as well as others (ref. 2), that the properties of pyrolytic, pressure- 
annealed graphite approach the properties of natural single crystals, and in discussing 
the band structure will assume the two are the same. 

The effective mass Hamiltonian describing the graphite energy bands was derived 
by Slonczewski and Weiss (SW) (ref. 5). The Hamiltonian # is 
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I 



0 


where 


and 


H 


13 


IT* 

H 13 


0 

E 2 

H 23 

- H 23 






(34) 

h ; 3 

TT* 

H 23 

E 3 

H 33 


_ H 13 

- H 23 

H 33 

E 3. 


e i= 

a + yj r + 

2 5 



E 2 = 

A - y^T + 

2 & 


(35) 



T = 2 cos 7 r£ (36) 

where £ is the coordinate along k z measured from the point K in the Brillouin zone 
as shown in figure 5. Figure 5 shows the Fermi surface which is located along the 
Brillouin zone edges HKH (ref. 2). Whether the charge carriers at A and B are elec- 
tron or hole (one must be electron, the other hole) is associated with the signs of the 
band parameters. These parameters are y Q , y 1 , y £> r 3 , y 4 , r 5 , A, and which 






Minority carrier 
pockets predicted 
by SW theory? 



Figure 5. - Brillouin zone of graphite. Fermi surfaces along 
H, K, H are labeled A and B. 
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can be determined from experiment. Magnetoreflection and deHaas - van Alphen experi- 
ments have given the best values so far (ref. 23). 


7 0 = 2. 85±0. 10 eV 
y x = 0. 31±0.02 eV 
r 2 = -0.0185 eV 
7 g = +0. 29±0. 02 eV 
r 4 = +0. 18±0. 02 eV 
r 5 = -0.0185 eV 
A = +0. 009 eV 
£ f = -0.022 eV 


-\ 



(37) 


The signs of these coefficients were determined by Schroeder et al. (ref. 6), and the 
constants with arrows change sign if the assignment of electrons and holes is reversed 
(see fig. 5). A direct and independent determination of signs is made from the results 
of the present experiments as will be discussed in the section "Consequences of Electron 
and Hole Carrier Assignment to the Energy Band Structure. " 

The factor distorts the Fermi surfaces from having circular cross section per- 
pendicular to the "C" axis. However, this has little effect o r the extremal cross sec- 
tional areas or effective masses. Consequently, y t wiU be ignored (set equal to zero) 

o 

in any later discussion involving the band parameters. 

In zero magnetic field the energy bands are described, neglecting y g, by 


S = - (e 4 + e 2 ) • 


| (Ej - E 3 ) 2 + 7q(1 - i/) 2 a 2 J 


1/2 


j. 

_4 


e = - (E 2 + Eg) ± - (E 2 - Eg) 2 + y 2 (l + v) 2 


1 


1/2 


(38) 


where E^, Eg, and Eg are given by equation (35). The radial coordinate, centered on 
the HKH edge, for directions perpendicular to k z is given by a, and v is given by 


^0 


(39) 
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where r is defined by equation (36). 

In figure 6 energy (from eqs. (35) and (38) is plotted as a function of the coordinate 
ij along the k direction for cr = 0. On the left in figure 6 the bands are plotted with 
the old signs (negative of those indicated by arrows in eq. (37)). On the right in figure 6 
the signs are those indicated in equation (37). Figure 6 also shows that £p is measured 
from the bottom (old notation) or top (new notation) of the conduction band at points H. 



Figure 6. - Energy bands for graphite against the £ coordinate along 
k (from eg. (35)). 


The relation between the band parameters and the extremal cross sectional area of 
the Fermi surface at the center of the zone edge (point K in fig. 5) is 


A K = <^2 - e F>< a - 2 ’l + 2 >5 - V 


(40) 


3a on> 


and the effective mass at the center of the edge is 


m K = ' — 


i±L(lA + , 1 + r 2 + r 5 -r F ) 


3 a o y o 


(41) 


With spin-orbit interaction included in the Hamiltonian, the tips of the ”A” Fermi 
surface sections shown in figure 5 break into a separate tiny pocket of carriers located 
at H as marked. The extremal cross sectional area A m of this surface, perpendicu- 
lar to the k_ direction is predicted (ref. 2) to be 

2D 
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A m = 


4ti£ f (S f - A) 

s 4*o 


(42) 


For a magnetic field present the energy levels were derived by McClure (ref. 24) 
and Inoue (ref. 25) from SW theory. For r 3 = 0 the Hamiltonian in a field has matrix 
elements given by 


JSf = 


E 1 

0 

H 13 

H i 4 

0 

E 2 

H 23 

H 24 

H 13 

H 23 

E 3 

0 

H 14 

H 24 

0 

E 3 


(43) 


where Ej, E 2 , and Eg are given by equation (35) and 




13 


\2 / 


H 


' 2 \ V2 (1 - ,)„ V 2 

^ / 

a/2 


14 


H 23 = ^ (! + y )( n + l ) l/2 

H 24 = ^ V 2 

2 = 3a o>l — 

° ° 2KC 

o 

and aQ is a lattice constant equal to 2. 456x10 centimeter. The secular equation 
from equation (43) is 


(44) 


"1 Q = 


E - Er 


E — E - E 2 


1(1 -V) 2 (l+l/) 2 . 


E - Er 


E - E^ E - E 2 


\_(l-v ) 2 (l + v) 2 _ 



(45) 
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EXPERIMENTAL 


Three different magnets were used for the present investigations: an 18. 5-tesla 
(1 T = 10 kG) liquid-neon- cooled steady state solenoid (ref. 26), a 10.5-tesla super- 
conducting solenoid, and a 4. 0-tesla transverse split pair superconducting magnet. The 
last magnet was used when measuring thermoelectric power as a function of angle be- 
tween the magnetic field and the graphite C axis. For the 10.5- and 4.0-tesla magnets, 
field measurements were accurate to within 0.2 percent. The samples were highly 
ordered pressure annealed pyrolytic graphite. The distribution of directions of "C" axes 
of the various crystallites was measured by X-ray scattering to be on the order of 2°. 
Transmission electron microscope studies found typical grain size to be about 10- 
micrometer diameter. The ratio of electrical resistance at 300 K to that at 4.2 Kwas 
1. 5 for sample PG 3, 9 for PG 4, and 9. 5 for PG 5. By comparison, natural single 
crystals have had ratios ranging from 12 to 37 (ref. 3). 

The sample geometry is shown in figure 7. For p^ and measurements, 

electric current flows along the sample length (±y direction) perpendicular to the field 



Figure 7. - Experimental geometry. Coordinate system defines components 
of tensors; z direction defined by field; y direction defined by sample 
length and current flow. 
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and voltages are recorded. The measured electric fields are E*, and from equation (3) 

"E* = pJ (46) 

when there are no temperature gradients. For adiabatic thermopower e^. and Nernst- 
Ettingshausen measurements, a constant heat flow w* is maintained in the sample 

and equation (1) shows 


E* = e’w* (47) 

To measure the thermal resistivity tensor components, equation (8) is used 

G = yw* (8) 

where G is the negative temperature gradient vector. To calculate the tensor elements 
of the isothermal thermopower e, we use equation (10) to obtain X from y and equa- 
tion (11) to obtain e from e\ 

Voltages are measured with a direct- current nanovoltmeter and recorded against 
field or angle on an x-y recorder. Temperature gradients are measured with an 
alternating-current Wheatstone bridge. Evanohm is used for heater wire to supply the 
constant heat current w*. The power from an Evanohm heater was shown to change by 
less than 2 percent in going from 0 to 18.5 tesla (ref. 27). 


EXPERIMENTAL RESULTS 
Magnetoresistance 

In figure 8 the magnetoresistance p and Hall resistivity p^ are plotted as a 
function of field to 10 tesla (1 T = 10 kG) at 1. 1±0.05 K. We find experimentally that 
py X « p^ and that quantum oscillations are obvious in both quantities. Using equa- 
tion (12) for fields parallel to the high symmetry C axis 


a 


yy 


yy 


p yy + 



(48) 


Thus maxima predicted by equation (13) for Oy. y result in minima for p yy when Landau 
levels coincide with the Fermi energy. Since 
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Figure 8. - Hall resistivity and magnetoresistance pyy as functions of 
field at 1. 1 K. Sample PG 5* scales are not the same. 
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Figure 9. - Oscillatory part of magneto resistance as function of field at 4. 2 
and 1.1 K. 
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where p means the oscillatory part of the resistivity. From equation (16) for the 
oscillatory conductivity 


p ocp° /5 sA 1/2 > Lnh^VTV 1 (2.M)- 1 / 2 (50) 

yy yy U F i / . to ^ kc o ) \ h 4 / 

L *M 

In figure 9 experimental results of the oscillatory part of the resistivity is plotted as a 
function of field to 4 tesla. At 4. 2 K the oscillations are undistorted up to nearly 2 tesla, 
and the oscillations are described by the M = 1 term in equation (50) for low fields. At 
1. 1 K the oscillations are spike shaped and the minima shifted to higher fields. It takes 
more than the M = 1 term in equation (60) to describe the oscillations. In the limit of 
near zero temperature and near perfect crystals, the resistance minima would be very 
sharp and drop nearly to zero as follows from equation (13). 


TABLE I. - COMPARISON OF HOLE MAJORITY CARRIER EXTREMA 
FOR LANDAU LEVEL CROSSING OF e F WITH PREDICTIONS 


OF SUGIHARA AND ONO (REF. 11) 


Landau 

Theoretical 

Thermopower 

Resistance 

Conductivity, 

level 

field position 

maxima 

minima 

°xy 


(Sugihara and Ono) 



minima 



Field, 1 



1 

3.6 

3. 66±0. 02 

3. 33±0. 02 

3. 64±0. 1 

2 

1.99 

1. 97±0. 03 

1. 88±0. 02 

1. 94±0. 04 

3 

1.42 

1. 39±0. 01 

1. 36±0. 02 

1. 39±0, 02 

4 

1.10 

1. 08±0. 02 

1. 01±0. 02 

1. 11±0. 05 

5 

— 

. 88±0. 02 

. 83±0. 02 

. 86±0. 03 

6 

— 

. 75±0. 03 

. 72±0. 03 

. 75±0. 03 

7 

— 

. 64±0. 03 

. 60±0. 03 

. 63±0. 02 
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TABLE II. - COMPARISON OF ELECTRON MAJORITY CARRIER EXTREMA 


FOR LANDAU LEVEL CROSSING OF WITH PREDICTIONS 


OF SUGIHARA AND ONO (REF. 11) 


Landau 

Theoretical 

Thermopower 

Resistance 

Conductivity, 

level 

field position 

minima 

minima 

°xy 


(Sugihara and Ono) 



maxima 



Field, 

T 


1 

6.4 

„ f7. 68±0. 05 

0 f7. 6*0. 1 

, f7. 4*0.1 



a 7. 0< 

a 7. 2< 

a 7.lJ 



(.6.92*0.05 

1.6. 8*0.1 

1.6. 8*0.1 

2 

2.95 

2.98*0.02 

3.0*0.05 

3.0*0.05 

3 

1.95 



1.86*0.02 

4 

1.54 



1.50*0.02 

5 

— 



1.22*0.02 


a Spin split Landau level. 


The field values for resistivity minima occurring at coincidence of S-p and Landau 
levels are predicted by Sugihara and Ono (ref. 11). These are listed in tables I and II for 
the quantum numbers n for the hole carriers and electron carriers, respectively. The 
resistivity minima from figure 9 for 1.1 K are also listed. As expected from the theory 
(IB) the resistance minima occur at slightly lower fields than Landau level crossings. 
Included in the table are the minima near 7 tesla shown in figure 8. At 1. 1 K the weak 
n = 2 minima due to the electron carriers is visible (figs. 8 and 9) at 3 tesla. The double 
minima near 7 tesla for the n = 1 electron carriers is double due to spin splitting of the 
Landau levels. This splitting will be discussed in the section ’’Spin Splitting. " The 
deHaas - van Alphen period of the minima in figure 9 is 0. 206±0. 005 tesla" h 


Hall Coefficient 

Figure 10 shows a plot of resistivity p , Hall resistivity p , and Hall conduc- 
tivity times H^, calculated from p,_ r and p._ using equation (12). The positions of 

yy y* 

extrema in all three quantities are marked by field value and Landau quantum number. 

As discussed previously Argyres (ref. 8) predicts that maximal extrema will result in 

<r T when Landau levels for electron carriers cross and minimal extrema will re- 
xy p r 

suit when hole levels cross £p. In figure 10, a X yH is plotted rather than a xy in 

order to best observe extrema, and has (spin split) maxima at 7 tesla, and maxima at 

2 

3.0 and 1. 8 tesla. This series of extrema is thus due to electron carriers, a H has 
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n = 3 n = 1 electron 



Figure 10. - Representative data for Hall resistivity, magneto-resistance 
and Hall conductivity as function of field at L 1 K. Sample PG 5. 


minima at 3. 64 and 1. 94 tesla, and this series is due to hole carriers. The field values 
correspond to crossings of Landau levels predicted by Sugihara and Ono (ref. 11) and 
these are listed in table I for holes, and table II for electrons. 

In figure 11 the measured p^ x is plotted as a function of field at 1. 1 K. A compli- 
cated series of superimposed maxima and minima is observed at low fields. A few of 
the extrema are marked by field and quantum number in figure 11. The best values from 
all data including figure 11 for these extrema are also listed in tables I and II. 

By comparing the ct extrema with resistance minima and Sugihara and Ono pre- 
xy 

dictions, we find that a minima belong to the deHaas - van Alphen series with period 

i 

0. 205±0. 005 tesla . Soule, McClure, and Smith (ref. 3) (SMS) have shown that this 



Figure II. - Representative data for Hall resistivity as function of field at 
1. 1 K showing low field detail in sample PG 5 and spin splitting near 
7 tesla. 
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deHaas-van Alphen period belongs to the Fermi surface section marked "A" in figure 5. 
This frequency assignment is also predicted by the SW band structure theory. Thus the 
Fermi surface sections at "A” in figure 5 are hole type carriers. 

The or maxima, when compared with Sugihara and Ono predictions (table n) are 
xy 1 

part of the deHaas - van Alphen series with period 0. 16±0. 01 tesla" and located at "K" in 
the Brillouin zone (fig. 5). The Fermi surface section B in figure 5 is thus the electron 
section. Sugihara and Ono theory and experiment are in accord for all data except for 
the n = 1 (spin split) electron level near 7 tesla. 

The major conclusions from the magnetoresistance and Hall effect results are: 

(1) The carriers marked A in figure 5 are holes and those marked B are electrons. 

(2) The Sugihara and Ono theory is within experimental uncertainty except for the 
n = 1 electron level which is found higher in field than predicted. 

(3) The Adams and Holstein theory for resistance in strong magnetic fields is quali- 
tatively correct in predicting (a) sharp maxima at low temperatures for o - , (b) maxima 
shift to higher field values as the temperature decreases, and (c) resistance minima 
(conductivity maxima) for both electron and hole carriers. 


Thermoelectric Power and Thermal Conductivity 

The adiabatic coefficient are experimentally measured, but the isothermal coeffi- 
cients from equation (1) with J = 0 must be calculated in order to compare with theory. 


E* 


= e’w* 


(51) 


so with w* in the y direction 


and 


E = e' w 

y yy y 


(52) 


E = €' w 
x xy y 


(53) 


(For relative coordinates see fig. 7. ) The magnetic field is maintained along the z di- 
rection. From equations (52) and (53), a measurement of Ey with constant Wy gives 

e ) the adiabatic thermopower, and E v gives e* the adiabatic Nernst-Ettingshausen 
yy x yx 

coefficient. To obtain the isothermal e tensor we use 


and 


e= e'X 

(ID 

A = f 1 

(54) 
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where y is the thermal resistivity tensor. For fields along the C axis of graphite 


y = 



r. 


xy 


r. 


yy 

o 



( 55 ) 


and 
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yy 


yy 2 2 

y yy + r xy 


(56) 


r. 


xy 


y* 2 2 

r yy + y xy 


(57) 


Experimentally we find y xy « y y for all ranges of magnetic field. Thus 


yy 


y. 


yy 


(58) 


and 




y. 


yy 


(59) 


The amplitude of Nernst-Ettingshausen coefficient oscillations were up to eight times 
larger than the thermopower. However, y xy was smaller than y yy divided by 50. 
Equation (10) shows that 

y = A _1 (10) 


where symmetry makes x^ = y^, y xy = -y^, and y xz = y yz - y zx = y zy = 0. Since 
y yy >: 50 y xy it follows from equations (10) and (11) that 


e 

yy 


e' x 
yy yy 


(60) 


and from equation (58) 
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(61) 


We have measured as a function of field to 10 tesla and find that it saturates to a 
constant value after increasing for about 0. 1 to 0. 2 tesla. This saturation is due to a 
freezing out of thermal conduction by electrons. As the field increases, conductivity by 
electrons decreases until it is unimportant compared with lattice conduction. Lattice 
conduction is independent of field. We did find very weak quantum oscillations in the 
thermal resistivity in the saturation region of magnetic field. These were in phase with 
the electrical resistivity oscillations, so were attributed to remnant electronic conduc- 
tivity and not to the phenomenon of phonon scattering from quantized electron states 
(ref. 28). 

Since is independent of field, equation (61) shows that 

£ yy (H > “ < 62) 


at a fixed temperature. Thus measurements of as a function of field are propor- 

tional to £yy as a function of field. In Ziman's (ref. 13) notation for Thermopower, 
Syy (see eq. (22)), 



(63) 


In figure 12 the thermoelectric power of sample PG 3 is plotted as a function of 
magnetic field to 2.4 tesla at 1.2 K. This figure shows that at low fields a very small 
frequency oscillation (minority carriers) is superimposed on the larger frequency (ma- 
jority carrier) oscillations. These oscillations could possibly result from carriers 
near the Brillouin zone boundary at points ”H" shown in figure 5. In figure 12 the 
minima in thermopower for the majority carrier oscillations are marked by arrows. 

In sample PG 3 (resistance ratio 1. 5) we studied the oscillation periods (P = 1/F). 
For the electron carriers the deHaas - van Alphen period is 0. 155±0. 005 tesla - * for 
field parallel to the C axis. These oscillation periods were studied over a wide range 
of angle between the field and the C axis and dominated at low field and low temperature. 
At 4. 2 K, oscillations from the hole Fermi surface dominated. Parallel to the C axis 
this period was 0.208±0.008 tesla . In figure 13 the electron and hole (majority car- 
rier) periods are plotted as a function of angle out to 70° to 80° from C. The solid lines 
are from the natural single crystal data of Soule, McClure, and Smith (ref. 3). The 
electron (lower curve) period found in PG 3 was slightly higher than the Soule, McClure, 
and Smith results but experimental errors overlap. The period against angle data of fig- 
ure 13 is characteristic of a nearly cylindrical Fermi surface, and Soule, McClure, and 
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Figure 12. - Thermopower S^. as function of field at 1. 1 K. Two carrier oscillations: majority electrons and 
minority carriers. Sample y ?G 3. 



Degrees between C axis and field 

Figure 13. - deHaas - van Alphen period P as function of angle between field 
and C axis. Majority carriers in sample PG 3. 
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Smith found the ratio of areas for the hole surface in natural single crystals to be 

A(90 ^ l = 12.1 (64) 

A(0°) 

where A(0°) is the extremal cross sectional area of the Fermi surface for a field paral- 
lel to the C -axis and A(90°) is for a field perpendicular to C. Within experimental er- 
ror we find that natural single crystals and pyrolytic graphite have identical areas for 
both electrons and holes out to ±70° from the C-axis. 

In figure 14 the deHaas - van Alphen frequency is plotted against the angle between 
the magnetic field and the C-axis for the same data as appear in figure 13. The fre- 
quency is proportional to the extremal cross-sectional area of the Fermi surface per- 
pendicular to the field. Figure 14 demonstrates how little change in area occurs for 
angles less than ±40° from C. From 60° to 90° the change in area is very rapid. 

In figure 15 the thermoelectric power of sample PG 4 (resistance ratio, 9) is plotted 
as a function of field to 10 tesla. The field values for electron and hole extrema are 
marked, as well as the quantum numbers associated with each peak. In the theory sec- 
tion, it was concluded that sharp positive peaks result from hole carrier Landau levels 
crossing £p. The series 3.63 tesla, 1.96 tesla, etc. , are thus due to hole Landau 
levels crossing £p. They have a deHaas - van Alphen period of 0. 206±0. 007 tesla - '*'. 

As discussed earlier, this period originates from the Fermi surface section marked A 
in figure 5. The A surface is thus a hole surface. The minima seen in figure 12 at 
7 tesla (spin split) are due to electron Landau levels crossing £p and show that the 
Fermi surface B in figure 5 is for electrons. This assignment of electrons and holes 
agrees with out Hall results. We have thus provided two independent and direct experi- 
ments to confirm the assignment of electrons and holes in the Brillouin zone. This as- 
signment agrees with the results of Schroeder, Dresselhaus, and Javan (ref. 7). 

In figure 16 the thermopower of sample PG 5 (resistance ratio, 9.5) is plotted as a 
function of field to 4 tesla. These data are the sharpest of the three samples. The 
same assignment of electrons and holes in the Brillouin zone results as with PG 4. An 
interesting part of these data is the line shape of the oscillations. The trace, as a func- 
tion of increased field, rises rapidly and plateaus coming down the high field side. The 
thermodynamic theory predicts the plateauing first and the sharp drop on the high field 
side as shown in figure 3. The same is true for the theories presented in references 18 
to 20. However, Horton (ref. 14) predicts the observed line shape, as shown in fig- 
ure 14. This lends support to the Horton theory. 

The observed spikes in thermopower occur at fields close to those predicted by 
Sugihara and Ono. The one exception is the spin split n = 1 electron level centered 
near 7 tesla. This peaks at a higher field than predicted. All of the observed crossings 


27 



deMaas - van Alphen frequency, T 



Figure 14. - deHass - van Alphen frequency as function of angle between field and C axis. Majority 
carriers in sample PG 3. 



Figure 15. - Thermopower against field at 1. 1 K for PG 4. Landau level 
crossings for both electrons (spin split) and hole carriers marked. 
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Figure 16. - Thermopower as function of field at 1. 1 K for sample PG 5. 


of Landau levels at are listed in tables I and II and compare favorably with Sugihara 
and Ono theory. 


Consequences of Electron and Hole Carrier Assignment to 
the Energy Band Structure 


Previous to the work of Schroeder, Dresselhaus, and Javan (ref. 7), the energy 
bands were as shown on the left in figure 6. The parameters yg, T5, and £jp were 
positive, and y ^ and A were negative. As shown on the right in figure 6, when elec- 
trons are located near the center of the zone edge (£ = 0), the coefficients yg an< ^ £p 
reverse sign. Our quantum Hall effect, and thermoelectric power results confirm that 
holes are at A and electrons are at B, as shown in figure 5. This makes y2 and 
negative. From equation (35) Eg then has negative values. The McClure-Inoue (refs. 
24 and 25) equation 


n + — )Q = 


E - Er 


E - Ej E - Eg 


Lu-i ') 2 d+d 2 J 


E - E , 


E - Ej E - E 2 


TV 


(l-u ) 2 (l+P ) 2 


Q 


2 \!/2 


1 


(45) 


can still be satisfied if the signs of A, y^, and y^ are reversed. From equation (35), 
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( 65 ) 


1 2 

E, = A + y 1 r + - 7 c F becomes -Er 


E« 


= a - y 1 r + ly 5 r 2 becomes -E^ 


= -y« r 2 


E 3 _ 2 7 2 J 


becomes -Er 


In equation (45) replace Ej by -Eg, Eg by -Ep E g by -Eg, and E by -E, and the 
equation is unchanged. Equation (40) for the cross sectional area predicts the same 
area in both assignments of sign. Equation (41) for the effective mass is also unchanged 
in value by the sign change. 


Spin Splitting 

Spin splitting of the n = 1 electron Landau levels are observed in the magneto- 
resistance (figs. 8 and 10), the Hall coefficient (figs. 8, 10, and 11), and the thermopower 
(fig. 15). For all of these, the separation of peaks is between 0. 6 and 0. 7 tesla. This 
is further confirmed by the observation of very nearly the same magnitude for splitting 
in three different pyrolytic graphite samples. In a natural single crystal, we have re- 
cently observed splittings in the n - 1 and n = 2 electron levels and in the n = 1 hole 
level. A study of consequences of the spin splitting is discussed in a separate paper 
(ref. 29). 


CONCLUSIONS 

The major findings of this study are 

1. Carriers marked A in figure 5 are holes and those marked B are electrons. 
This assignment is made from Hall effect and independently from the thermopower 
quantum resonances. 

2. The theory of Sugihara and Ono is within experimental uncertainty of predicting 
the observed positions of Landau level crossings of ip. The one exception is the n - 1 
electron level, which is found higher in field than predicted. 

3. The Adams and Holstein theory for resistance of a degenerate conductor in a 
strong magnetic field, is correct in predicting sharp maximas for at coincidence 
of a Landau level and the Fermi energy for both electron and hole type carriers. Shifts 
of peaks to lower fields at higher temperatures is observed. 

4. A thoery by Horton appears to predict the correct line shape of quantum oscilla- 
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tions in the thermopower at high fields. A Dingle temperature of 1.5 K gives a good fit 
between Horton theory and experiment. Other theories were found not to satisfactorily 
explain experimentally observed wave shapes. 

5. We have made the first observation of spin splitting of Landau levels in graphite. 
The n = 1 electron Landau level is spin split by 0. 6 to 0. 7 tesla at 7 tesla. This pro- 
vides a measure of the spin orbit splitting at the point K in the center of the Brillouin 
zone edge. 

6. The thermal resistivity saturates to a constant value independent of magnetic 
field strength for fields above about 0. 2 tesla. 

7. Pyrolytic and natural single crystals produce the same deHaas - van Alphen fre- 
quencies as a function of angle between the field and the C axis. This lends support to 
the use of the single crystal band theory for pyrolytic samples. 

8. A low deHaas - van Alphen frequency (minority carrier) oscillation is observed 
which could possibly originate from carriers in the points H in the Brillouin zone. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, June 18, 1970, 
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